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ABSTRACT 

About two-thirds of present-day, large galaxies are spirals such as the Milky Way or Andromeda, 
but the way their thin rotating disks formed remains uncertain. Observations have revealed that half 
of their progenitors, six billion years ago, had peculiar morphologies and/or kinematics, which exclude 
them from the Hubble sequence. Major mergers, i.e., fusions between galaxies of similar mass, are 
found to be the likeliest driver for such strong peculiarities. However, thin disks are fragile and easily 
destroyed by such violent collisions, which creates a critical tension between the observed fraction 
of thin disks and their survival within the A-CDM paradigm. Here we show that the observed high 
occurrence of mergers amongst their progenitors is only apparent and is resolved when using morpho- 
kinematic observations which are sensitive to all the phases of the merging process. This provides 
an original way of narrowing down observational estimates of the galaxy merger rate and leads to a 
perfect match with predictions by state-of-the-art A-CDM semi-empirical models with no particular 
fine-tuning needed. These results imply that half of local thin disks do not survive but are actually 
rebuilt after a gas-rich major merger occurring in the past nine billion years, i.e., two-thirds of the 
lifetime of the Universe. This emphasizes the need to study how thin disks can form in halos with 
a more active merger history than previously considered, and to investigate what is the origin of the 
gas reservoir from which local disks would reform. 

Subject headings: galaxies: evolution, galaxies: formation, galaxies: interactions, galaxies: kinematics 
and dynamics, galaxies: general 



1. INTRODUCTION 

Over the past years, observations and theory have been 
progressing significantly in understanding how the first 
large disk galaxies could form. Near-infrared (rest-frame 
optical) imaging from the HST reveals that, while ~40% 
of z^l-3 most massive (i.e., M ste iiar > 5 x 10 10 M Q ) 
galaxies are relatively compact and quiescent, ~60% 
of them are actively for ming stars and sho w quite 
extended stellar structur e s dWeinzirl et al.ll2011t see als o 
Ivan Dokkum et all 120111 : IForster Schreiber et al.l 120111) . 
In a la rger range of mass (i.e., ~lO 9 -lO 11 M0L lLaw et al.l 
(|2012l ) suggested that the stellar phase of z~1.5-3.6 
galaxies is probably distributed as triaxial ellipsoids 
rather than thick disks. Integral field spectroscopy now 
routinely provides us with spatially-resolved kinematics 
of th e sub-population of z~ l -3 em is sion line galax- 
ies (IForster Schreiber et al.l 120061: IBournaud et al 



20081: Ivan Starkenburg et all 12001 iWright et al 



200a iForster-Schreiber et all 120091: lLaw et al l 1 20091: 



Epinat et all 2009t iLemoine-Busserolle et al.l bold 



Lemoine-Busserolle fc L amarcillel 120101 : iContini et all 
2012t lEpinat et al.l I2012D . Actively forming galaxies 



at these epochs reveal a relatively balanced mix of 
obvious (major) mergers, relatively compact dispersion- 
supported objects, and rotating extended structures, 
all being much more turbulent in their gaseous phase 



comp ared to local galaxies (jForster-Schreiber et al.l 
2009). On the theoretical side, cosmological simula- 
tions have predicted that streams of cold gas would 
occur in the halo s in which these galaxies are ex- 
pected to inh abit (iDekel et al.l 120091: iKeres etldl 120091: 
l Agertz et "all 120091: [Brooks et all 120091: iCeverino et al.l 
120101 : Ivan de Voort et all 1 2 1 lh . Such streams would 
be responsible for bringing a continuous amount of 
high angular momentu m gas in these halos over a rela- 
tively short timescale ()Kimm et al.l 120111 : iPichon et al.l 
1 2 1 If ) . Because the amount of mass brought that 
way is expected to be considerable, this could in 
principle pressurize the central forming di s k and 
result in clump fragmentation (iDekel et al.1 [2009b; 
IBournaud et~aTl 120091: iGenzel et al .l 1201 lft. as ubiq- 
uitously ob served at these epochs (|Elmegreen et all 
120051 I2009HH ). Such "clumpy galaxies" are now more 
and more considered as an important step of galaxy 
formation. However, it remains unclear whether the 
lifetime of such clumps is long enough to allow them 
to survive fee dback effects and coa l esce to form a 
centra l bulge (Elmegr een et al.1 120081: IBournaud et all 
2009 b: iMurrav et all 120101: iKrumholz fc Dekell 120101: 
Genel et all 120101: (Forster Schreiber et all 1201 lbi 
IPowell et alll2011l:lHopkins et alll201lD . It remains also 
unclear how thin disks can form from the high angular 
momentum gas which is expected to be accreted into 
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the halos. Finally, the smoking gun detection of cold 
gas accretion is still missing, alt hough different options 
are being ex tensiv ely discussed (Diikstra fe Loebl 120091: 
20101 iGoerdt et al 1120101 iBenson fe Bowerl 
1 120111: iFumagalli et al l 
Le Tiran etall 120111 : 



Steidel et al. _ 

2011 1 : iFaucher-Giguere fe Kere 



201lt iKimm et all l20lTbT 



Giavalisco et al.ll2011l: [Stewart et al.ll2011llbf ). 



In contrast, it remains more uncertain how these z~2 
thick rotating disks evolved down to z~0, and how the 
local thin disks formed. Indeed, the Hubble sequence of 
galaxies such as the Milky Way (M steUar = 10 10_11 M Q ) 
reveals a predominance (72%) of rotationally-supported 
disks in the present-day Universe ([Delgado-Scr rano et al.l 
2010). Using the Cosmological Principle, their progen- 
itors at z~0.6, i.e., six billion years ago, can be se- 
lected and their morphology studied using the Advanced 
Camer a for Survey on-board the Hubble Space Tele- 
scope (Dclgado-Se rrano et al.l [2010h . Using methodolo- 
gies and rest-frame luminosities that mimic those used 
for local galaxies, the likeliest an cestor of the Hubble se- 
quen ce can then be established dDelgado-Serrano et al.l 
2010). Half the progenitors of spiral galaxies show 
peculiar morphologies that ex c lude them from this 
seque nce (|van den Berghl 120021 : iDelgado-Serrano et al.l 
I2010D . Most of them are starbursts, whose emer- 
gence at z~0.6 refle c ts that of emission line galaxies 
(|Hammer et"aH Il997t iMignoli et all 120051) . As part of 
the IMAGES survey, the large-scale spatially-resolved 
kinematics of a representative subsample of 63 emis- 
sion line galaxies at z=[0.4-0.75] was obt ained using 
the GIRAFFE multi-IFU spectrograph (jFlores et all 
[200l iPuech et al.l [20061: lYang et al.l [2001 . Most pe- 
culiar galaxies also reveal peculiar large-scale motions 
(jNeichel et al.l [20081 ). Comparing the morpho- kinematic 
properties o f the subsample o f the 33 galaxies lying into 
the CDFS (|Yang et al.l [2001 ) with those of a grid of 
simple major merger models, convincin g matches were 
found in about two-thirds of the cases (|Hammer et al.l 
2009). This leads to at least one third of the progeni- 
tors of present-day spirals potentially involved in a ma- 
jor merger. Sin ce major mergers can eas i ly destroy thin 
rotating disks (iToomre fe Toomrel 119721: iToomrd Il977t 
Barnesj|1988l 119921: iHernquistj 119931: iBarnesfe Hernquistl 
19961: iToth fc Ostrikerl 1199$ Hopkins et al.l 120081) . this 



creates a critical tension between the large fraction of 
present-da y disks and their survival within the A- CDM 
paradigm ([Stewart et~aTll2008l : IWeinzirl et al.ll2009f ). 

Further investigating the "disk survival issue" requires 
an accurate estimation of the evolution of the merger 
rate, i.e., the fraction of galaxy involved in mergers per 
unit time. Several techniques were developed over the 
past years to try to quantify accurately the evolution 
of the merger rate, such as the pair techn i que (e.g., 
Le Fevre etall 120001 iKartaltepe et~aT1 120071: iLin et al 



2008; Rawat et al. 2008; Ryan et al. 2008; Bundv et al 



2009? Ide Ravel et alj |2009l. 12011 " or the morpholog- 
ical technique, either through non- parametric (e.g., 
Conselice et all 120031 : lLotz et all 120081: iConsel i ce et all 
2009L iShi et al.l 120091: lLopez-Saniuan et al.l l2009D or vi- 
sual (e.g.. Uogee et al.l 120091: iBridge et all I2010D clas- 
sifications. Alternatively, the amplitude of the two- 
point correla tion function can also be used to in fer a 
merger rate (|Bell et all [20061: IRobaina et al.ll2010h . All 
these techniques are very dependent on the visibility 



timescale of the specific merging phase sampled by the 
observations, which requires to be estimated using ad- 
ditiona l numerical simul ations and/or analytical argu- 
ments (jLotz et al.ll2008bl ). Accounting for the exact dis- 
tributions of merger mass rati os or orb its involved is also 
challenging (jLotz et al.ll2010l lb1. 120111 ). As a result, ob- 
servati onal estimates spread over almost on e decade at 
z-0-1 (jHopkins et alll2009l; lLotz et alll20ll . 

Recent progresses combine numerical simulations with 
transfer radiation to derive a weighted average visibility 
timescale for each observational technique, using predic- 
tions for the distribution of mass ratios and ot her merger 
param eters from galaxy formation models (|Lotz et al.l 
1201 ID . Taking these aspects into account, the discrep- 
ancy between the different merger rates derived in the 
literature is significantly reduced, while remaining differ- 
ences are attributed t o differences in ra nge of mass ratio 
and sample selection ([Lotz et al.l 120111 ). This illustrates 
how crucial it is to carefully accounts for these aforemen- 
tioned effects. The present study puts this effort one step 
forward: we show that by combining spatially-resolved 
morphological and kinematic data as those provided by 
the IMAGES survey, one can date back the epoch where 
observed galaxy mergers were still in pair, and then de- 
rive the evolution of the fraction of mergers out to z~1.5, 
in a consistent way as a function of redshift. Such data 
also allow us to estimate relevant visibility timescales 
that are by construction consistent with the distribution 
of merger mass ratios and orbits of the galaxies in the ob- 
served sample. The derived evolution of the merger rate 
out to z~1.5 is found to be within a factor 2-3 of theo- 
retical predictions from a state-of-the-art A-CDM semi- 
empirical model. 

This paper is organized as follows. Sect. 2 presents the 
data used in this paper, which come from the IMAGES 
survey of z~0.6 emission line intermediate-mass galaxies. 
In Sect. 3, we argue why large morpho-kinematic pertur- 
bances most likely result from major mergers. Sect. 4 
presents the major merger fractions and rates derived 
from these data. A comparison with state-of-the-art 
semi-empirical model of galaxy formation is presented 
in Sect. 5. Finally, a discussion is given in Sect. 6 and 
conclusions are drawn in Sect. 7. Throughout the paper, 
magnitudes are given in the AB system, an d a WMAP- 
7 cosmology is used (jKomatsu et al.l [201 ID . with (fl m , 
fl A , Q b , h, cr 8 ) = (0.275, 0.725 0.0458, 0.702, .816). We 
adopt the diet Salpeter IMF ([Bell et all 120031 ) through- 
out this paper. Note that all IMF-dependant quantities 
cited from other studies were converted to this IMF. 

2. MORPHO-KINEMATIC DATA OF Z=[0.4, 0.75] 
GALAXIES: SUMMARY OF THE IMAGES SURVEY 

2.1. The representative IMAGES- CDFS sample 

In the frame of the IMAGES survey, a sample of 35 
emission line galaxies at z=[0.4-0.75], i.e., about six bil- 
lion years ago (mean/median redshift of z=0.6/0.65), 
were observed in t he Chandra Deep Field South (CDFS) 
(|Yang et al.l [20081) . This subsample of IMAGES is of 
particular interest sinc e it benefits from homogeneous 
imaging and 3D data (|Yang et all 120081 : iNeichel et all 
120081 ). IMAGES targets were selected to have I AB < 
23.5 and EWq{[OII\) > 15 A to guarantee their detection 
by th e multi-IFU spectrograph FLAMES /GIRA FFE at 
VLT (jRavikumar et al.|[2007t TFlores et alll2006l) . These 
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galaxies were selected to be intermediate-mass galaxies 
(Msteiiar = 10 10 - lO n M using a diet Salpeter IMF 
and simplified prescriptions to convert J- ban d luminos- 
ity an d B-V colors into stellar mass, following iBell et all 
2003). This range of mass is of particular interest since 
mos t of the star formation at z<l occurred i n such galax- 
ies jHammer et all 120051: IBell et all 120051: IZheng et all 
1200711 Since the fr action of E/SO did not rais e be- 
tween z~0.6 and z=0 (jDelgado- Serrano et al.|[20To() . z<l 
intermediate-mass galaxies are the likeliest progenitors of 
local spirals (see also App. B). 

From the full IMAGES sample, one galaxy was re- 
j ected because it turned out not to be a starburst galaxy 
(|Hammer et al. 2009; Y ang et aT]|2009f) , and another one 
because the HST/ACS images were corrupted and no 
morphological analysis could be conducted. The re- 
sulting sample is fully represent ative of the J-ban d lu- 
minosity function at z=0.5-1.0 (|Yang et al.l 12008). At 
these redshifts, the near-infrared luminosity functions of 
the global population an d of the blue sub-pop ulation of 
galaxies are very similar (|Cirasuolo et al.ll2007ft . At first 
order, emis sion-line and blue galaxies s ample the same 
population (Delga do-Serrano et a l. 2010), which in prin- 
ciple guarantees that the IMAGES-CDFS sample is rep- 
resentative of the star forming population of galaxies at 
these epochs. 

The total SFR of all the 33 remaining galaxies in the 
IMAGES-CDFS sample was estimated as the sum of the 
unobscured SFR measured from the rest-frame UV lumi- 
nosity at 2800A combined to the dust-reprocessed con- 
tribution, which was estimated from the I R luminosity 
as me asured using 24 ^m Spitzer fluxes (see lPuech et al.l 
[20101 for details). We find a mean SFR of ~12M /yr, 
which compares well with other UV+IR measurements 
in larger s amples at a simi lar range of mass and red- 
shift (e.g.. IJogee et alj |2009). In Fig. [TJ the location of 
the IMAGES-CDFS galaxies in the SFR-M steHar plan is 
compared to the sequence of starbursting galaxies ob- 
served by iNoeske et all (|2007ft at z=0.45-0.7. As ex- 
pected, the IMAGES-CDFS sample comprises Luminous 
InfraRed Galaxies, which populates th e high-SFR tail 
of th e distribution at SFR>14M Q /yr (|Le Floc'h et al.l 
2005), and less obscured (UV) starbursts, which pop- 
ulate the low-SFR tail of the distribution. This figure 
reveals that most (91%) IMAGES-CDFS galaxies lie well 
within the observed SFR-M ste z; ar . sequence at these red- 
shifts. The star formation rate density in the IMAGES- 
CDFS sample is found to be log {p s fr) ~-l-4 (where p s f r 
is expressed in Af Q .j/r _1 .Mpc -3 ), also in good agree- 
ment with expectati ons from larger su rveys that include 
IR data at 24 /jm (|Zheng et all 12007ft PI This confirms 
that the IMAGES-CDFS sample is representative of star- 
forming galaxies at z~0.6. 

1 Apparent discrepancies of up to a factor two can appear in 
comparison with other studies based on different survey areas, stel- 
lar mass cuts, and methods for estimating the SFR. For instance, 
IJogee et al.l 12003) found log (p s f r ) =-1.4 at z~0.6 for galaxies 
with M-stellar > 1O 9 ' 1 M0 and log(p s f r ) =-1.9 for galaxies with 
MsteMar > 1 10 5 Mp) (using a d iet Salpeter IMF). Using a larger 
area survey, IZheng e t al. ( 20Q3) found that for the mass range 
M s t e (j ar = 10 1015 — 1O 1115 M0, i.e., very close to the range of mass 
considered in the present study, log (p a f r ) ranges from -1.56±{5;{5g 
at z=0.4-0.6 to -1.34±0.05 at z=0.5-0.8, which is consistent with 
the value found in the IMAGES-CDFS sample. 




10.5 11.0 
Log(M stell0f /M 

sun/ 

Figure 1. The SFR-Stellar mass sequence in the IMAGES-CDFS 
sample (black squares). SFRs and stellar masses were rescaled to 
a diet Salpeter IMF. The grey region represents t he z=0.45-0.7 
SFR-S tellar mass sequ ence and its ±l-o sc atter from lNoeske et all 
12007ft , as reported bv lDutton et al.l 12OT0T) . 

Gas masses were estimated from t he total SFR by in - 
verting the Schmidt-Kennicutt law (|Puech et al.ll2010l) . 
All galaxies were found in a range of total baryonic mass 
log (M b /M Q ) =10.2 to 11, or log (M steUar /M Q ) in the 
range 10 to 11. The stellar-mass density in the IMAGES- 
CDFS sample is found to be log {p s teiiar) ~8.1 (where 
patellar is expressed in M & .Mpc~ 3 ), which is in relatively 
good agreement with expectatio ns fr om larger survey s 
of blue galaxies (jArnouts et all 120071; [Bell et all 12007ft . 
within expected uncertainties. 

2.2. Morpho-kinematic classification 

Detailed color morphology (at a ^500 pc resolution 
scale) and large-scale kinematics (at a ~7 kpc resolu- 
tion scale) were obtained thanks to the Hubble Space 
Telescope, and 3D spectroscopy using the GIRAFFE 
multi-IFU spectrograph, respectively. Galaxies were 
classified independe ntly according to their morpholog y 
(|Neichel et alJ 12008ft and kinematics (|Yang et all 12008ft . 
based on similarities with the properties of local galaxies. 

The kinematical maps allowed us to classify the 
IMAGES- CDFS galaxies in to three distinct kinematic 
classes (see lYang et al.ll2008l ). Rotating disks [RDs] were 
those for which the Velocity Field [VF] showed an order 
gradient, a dynamical axis aligned with the morpholog- 
ical axis, and a velocity dispersion peak located at the 
dynamical center. For these objects, it was required that 
the sigma peak location and amplitude can be repro- 
duced by modeling the velocity dispersion map from the 
velocity gradient observed in the VF. At the GIRAFFE 
spatial resolution, it is indeed expected that most of the 
rotation curve gradient falls into on GIRAFFE pixel, 
which results in an velocity dispersion peak at the dy- 
namical center of the object. Conversely, galaxies with 
complex kinematics [CKs] were those for which no or- 
dered gradient was observed in the VF, or when a sig- 
nificant misalignment between the kinematical rotation 
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axis and the morphological principle axis was detected. 
Perturbed rotators [PRs] gathered objects with an order 
velocity gradient in the VF, but for which the velocity 
dispersion peak (location and/or amplitude) could not be 
reproduced from the VF, which indicates that this peak 
cannot be attributed to rotation. This classification was 
tested against quantified criteria and considering all un- 
certainties associated to kinematic measurements, which 
makes it objective a n d rep roducible, as described exten- 
sively in lYang et all (j2008[ ). 

The detailed color morphology of all the 33 galaxies in 
the IM AGES-CDFS sample was studied bv lNeichel et al.l 
(2008). Using a reproducible decision tree relying on 
parametric parameters and visual inspection of images 
and color maps, galaxy morphologies were classified into 
six different classes: Compact systems, Mergers, Pecu- 
liar/ Possible Mergers, Peculiar Tadpoles, Peculiar Ir- 
regulars, and Spiral disks. In particular, it was required 
for a rotating disks to show redder colors towards the 
center as observed in local galaxies. This classification 
was done by three independent classifiers using a decision 
tree, which limits the subjectivity and makes it repro- 
ducible. The comparison between the morphological and 
kinematical classes show a very good agreement, with 
more than 80% of spiral disks having rotating disks, and 
more than 90% of peculiar galaxies having complex or 
perturbed kinematics. 

Since perturbations in the (ionized) gaseous phase gen- 
erally correspond to perturbations in the stellar phase, 
the result of both the morphological and kinematical 
classifications can be ro bustly synthesized int o a morpho- 
kinematic classification (lHammer "et~aLll2009h . With this 
classification, rotating spiral disks are those classified 
both as rotating disks from the kinematic point of view, 
and as spiral disks from the morphological point of view. 
Non-relaxed systems are those that show a complex kine- 
matics and a peculiar morphology. Semi-relaxed systems 
are galaxies that show either some rotation from their 
kinematics (being classified as either a rotating disk or 
either a perturbed rotation) but with a peculiar mor- 
phology, or galaxies that do not show rotation in their 
kinematics (being classified as having a complex kine- 
matics) but with a spiral disk morphology According 
to this classification, 80% of z~0.6 sta rburst galaxies 
are n on-relaxed or semi-relaxed systems (jHammer et al.l 
2009), which corresponds to ~50% of all intermediate- 
mass, z~0.6 galaxies, since 60% of inte rmediate-mass , 
z~0.6 galaxies are inde ed starbursts (Ha mmer et al.l 
[19971 IMignoli et al"][2005h . The HST/ACS images and 
kinematic maps for a representative sub-sample of the 
IMAGES-CDFS can be found in Fig. [6] (see App. A). 

3. ASSOCIATING KINEMATIC DISTURBANCES WITH 
PHYSICAL PROCESSES 

Stellar kinematic disturbances are indubitable evi- 
dence for an out-of-equilibrium dynamical state. In high- 
redshift galaxies, deriving spatially-resolved stellar kine- 
matics remains out of reach, and one has to rely on gas 
kinematics as a best proxy. This implicitly assumes that 
gas kinematics traces the gravitational potential, and 
that no other mechanisms are perturbing gaseous mo- 
tions, such as, e.g., outflows. In the IMAGES sample, 
large-scale o utflows have been detected in only a hand- 
ful of cases (|Puech et al.ll2010t ). Using FORS2 slit spec- 



troscopy, iRodrigues et al.l (|2012l ) further estimated that 
only 8% of intermediate-mass galaxies at z^0.6 have sig- 
nificant outflows. This suggests that gas kinematics is 
indeed a reliable tracer of the dynamical state of galax- 
ies, at least on large spatial scales. Several dynamical 
processes are known to be a source of kinematic distur- 
bances. Some of them are local in space, while others 
are global. Relative to the equilibrium state of rotation 
(as expected in the intermediate mass range since more 
massive galaxies arc expected to reach a different equilib- 
rium state supported by dispersion), local disturbances 
are either due to minor mergers or internal instabilities 
such as clump formation, while global disturbances can 
be due to major me rgers or cosmological gas accretion 
(jHopkins et ai'll2010h . 

3.1. Local disturbances 

In the IMAGES sample, kinematic perturbations are 
detected over a >7 kpc spatial scale, which corresponds 
to the GIRAFFE-IFU spatial resolution element (i.e., 
two spaxels). They generally affect the galaxies across 
a significant fraction of their diameter, which, in princi- 
ple, excludes de facto any mechanism resulting in local 
kinematic disturbances. In ord er to quantify this, we de- 
rived the Jeans scale length Lj (|Bournaud et a l. 2009) in 
the IMAGES-CDFS, which is found to be 0.7±0.1 kpc. 
Lj increases to 1.0±0.4 kpc i n the sub-sam ple of clumpy 
galaxies that show rotation (jPuechl l2010( ). Thus, such 
Jeans instabilities should remain largely undetected at 
the GIRAFFE IFU spatial resolution. In the bottom of 
Fig. [5]one of the most clumpiest galaxy in the IMAGES- 
CDFS sample. This galaxy shows at least seven clumps 
spread over its disks and is classified as CK. However, the 
kinematic perturbations are mostly large-scale, with a 
misalignment between the morphological and kinematic 
axis as defined by the location of the maximal and min- 
imal rotation velocities, and with a velocity dispersion 
peak that cannot be r eproduced b y the r o tation model 
inferre d from the VF I Yang et al.l (|2008D ; iPuech et all 
(2010). This confirms that even in the most disturbed 
clumpy galaxies, the kinematic perturbations cannot be 
associated with the position of the clumps, i.e., they are 
not local but global. In addition, clump instabilities are 
limited to a fraction of 20% of intermediat e- mass galax - 
ies at z^0.6 that are effectively clumpy (|Puechl [2010). 
which also rules out clump instabilities as a driver for 
the strongest kinematic perturbations. 

Regarding minor mergers, such processes are also ex- 
pected to lead to local kinematic perturbations only. An 
obvious case of minor m erger with a mass ra tio of ~1:18 
was studied in detail by IPuech et al.l (|2007[ ): only a lo- 
cal signature is found in the velocity dispersion map (see 
top panel in Fig. [2]). Even if we cannot exclude that 
in very specific cases such local processes could result 
in large kinematic disturbances such as those observed 
in IMAGES galaxies, we can safely conclude that they 
cannot be the main driver for the majority of galaxies, 
at least for the CK galaxies. This does not mean that 
minor mergers are not ongoing in the IMAGES sample, 
but given the spatial resolution of GIRAFFE-IFU data, 
they probably remain mostly undetected. 

3.2. Global disturbances 
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Figure 2. Examples of large-scale kinematic maps obtained in the frame of the IMAGES-CDFS survey with FLAMES/GIRAFFE at 
VLT. From left to right: HST/ACS z-band image superimposed with the GIRAFFE IFU grid (0.52 arcsec/spaxel), velocity field, velocity 
dispersion map, 5x5 linearly interpolated velocity field, and 5x5 linear ly interpolated vel ocity dispersion map. Top panel: J033226.23- 
274222.8, which is an example of a minor merger detected at z~0.67 (see[Pucch ct al. 2007). Bottom panel: J033230. 57-274518. 2, which is 
an example of a z~0.69 clumpy galaxy with Complex Kinematics. The local perturbations due to minor mergers or local instability remain 
undetected in the velocity field at the GIRAFFE spatial resolution (see text). 



Regarding the impact of processes resulting in global 
disturbances, it is instructive to consider the range of the 
relative variations of the gravitational potential (com- 
pared to the dynamical equilibrium state) A (/>/</>, in 
which galaxies can be considered as isolated and sub- 
ject to secular evolution, as a function of mass and 
redshift ([Hopkins et all 120101 ) . We estimated A</>/0 for 
galaxies in the I MAGES-CDFS sam ple using rotation 
kinematic models ([Puech et al.l 120081 ). For this, we de- 
rived the pixel-by-pixel relative difference between the 
observed VF and each galaxy model, and considered 
the sum AV/V of this map over all pixels, weighted by 
the observed S/N in each pixel. The quantity 2. AV/V 
can then be used as a proxy for A(j>/<f). We also de- 
rived the dynamical timescale td yn f° r each galaxy in 
the IMAGES-CDFS sample using as a proxy the ratio 
between the gas radius and the rotation velocity as de- 
termined bv lPuech et al.l (|2010f ). In galaxies showing ro- 
tation (RDs+PRs), we found that A<£/<£=1.2±0.4 (me- 
dian value) and ^j,n=60±3 Myr, while for CK galax- 
ies, A</>/(/>=3.1±1.0 and td !/ „=60±24 Myr. Uncertainties 
correspond to 1-er error-bar, and were estimated using 
bootstrap resampling. 

In this regime (i.e., tdyn <]10 Myr and Acj)/(j> >1), the 
two dominant physical processes expected to drive grav- 
itational perturbations in this range of mass and red- 
shifts are secular evolution (i.e., instabilities occurring 
in a galaxy) and the direct dynamical evolution of in- 
dividual substructures due to some previously initiated 
perturbation, which can ev entually dominate th e evolu- 
tion of the whole system (jHopkins et al.l [20101) . Large 
global perturbations (i.e., A<fi/<f> >1) are typically trig- 
gered by major mergers. The characteristic timescale for 
secular evolution can be estimated using the circular disk 
period i.e., Pd — 2ttR/V c , with R the disk radius and V c 
the circular velocity. We used ag ain gas radiu s and the 
rotation velocity as determined bv lPuech et al.l (|2010f ) as 
proxies for these two quantities. The timescale for dy- 
namical friction associated with the dynamical evolution 
of substructure was estimated following the simplified 
prescription given by [Hopkins e t al.l ([2010D : 



tdf 



0.2 



S Q ln(A) 



where So = A<fi/<fi, and A=l+l/5n. This timescale quan- 
tifies how fast the dynamical evolution of substructures 
evolve on their own, once they have been triggered. For 
instance, a merging galaxy will trigger instabilities in the 
main progenitors that will subseq uently evolve on thei r 
own, on a timescale given by tdf (jHopkins "eTaT][20Trll) . 
The dominant process at given amplitude of perturbation 
is the one associated with the shortest timescale. 

In Fig. [3[ we plotted the expected range of value of 
these two regimes at z~0.6 in the IMAGES-CDFS sam- 
ple. Secular evolution is found to be dominant at small 
relative variations of the gravitational potential A(f>/4>, 
while the influence of (major) mergers dominate above 
A<j)/4> >0.3. More specifically, isolated galaxies should 
be found in the upper left region of Fig. [3l close to the 
grey region (as well as large local disks such as M31), 
while on-going mergers and mergers remnants should be 
found at Acj)/(j> ~1 and close to the red region. Com- 
paring with these theoretical expectations, CK galax- 
ies lie closer to the individual dynamical perturbation 
regime (i.e., the merger-driven regime), which is consis- 
tent with (major) merger-driven kinematic disturbances 
in these galaxies. Galaxies showing rotation also lie in 
this regime, although with a smaller range of gravita- 
tional perturbation, which is consistent with these galax- 
ies being more (temporally) advanced merger remnants, 
as sugge sted by their other morpho- kinematic proper- 
ties (see lHammer et al. 2009| and Sect. 4.1). On over- 
all, z^0.6 rotators (RDs+PRs) lie closer to the expected 
merger-driven regime, which probably reflects the fact 
that they are closer to dynamical equilibrium than CK 
galaxies, in which gravitational perturbations are about 
three times larger. This suggests that any other known 
process (gas accretion or secular evolution) is very un- 
likely to be responsible for the high level of kinematic 
disturbance observed on large spatial scales, at least for 
the majority of the IMAGES-CDFS sample. 

3.3. Comparison with major merger models 

In this context, it was explored whether hydrodynam- 
ical major merger models for the 27 galaxies in the 
IMAGES-CDFS sample which are not relaxed morphog- 
ically and/or kinem atically could genera lly account for 
their perturbations (|Hammer et al.ll2009U bl: iPuech et al.l 
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Figure 3. Dynamical timescales (in Myr) for the evolution of sec- 
ular or self-amplified perturbations as a function of the amplitude 
of the gravitational perturbation. The grey region represents the 
expected 2-cr uncertainty range around the median timescale for 
secular perturbations, while the red region represents the expected 
2-cr uncertainty range around the median timescale for individual 
perturbations to evo lve on its own (as approximated by Eq. (5) of 
IHopkin s ct al. 20101. see text). The blue region is the 2-cr uncer- 
tainty range around the median value obtained for rotating objects 
at z^0.6, while the brown region represents the same for objects 
with a complex kinematics. 



[2001 iPeirani et all [20091: iFuentes-Carrera et al.l [20Toh . 
Morphological and kinematic maps from these simula- 
tions were produced and degraded to the spatial res- 
olution of the observations. Such maps were then ex- 
plored visually as a function of time and viewing angles 
to determine the best match to observations. The so- 
determined best model was graded by three different ex- 
aminers (independent of the one w ho determined the bes t 
model) using specific constraints (jHammer et al.l 12009). 
Secure major merger candidates (18/27) were identified 
as galaxies having a model that fitted observations with 
a good level of confidence (i.e., grades > 4/6). Exam- 
ples of such models are shown in App. A. These candi- 
dates can be used as a basis to estimate the fraction of 
mergers and the merger rate at different epochs as de- 
tailed in Sect. 4.1. Interestingly, the modeling exercise 
lead to a median merger mass ratio /i=3.03±0.34 in the 
IMAGES-CDFS sample, while theoretical expectations 
from variations of the gravitational potential give consis- 
tent predictions with A(j>/<f> ~ jU=3.1±1. in CK galax- 
ies (se e above). We refer the reader to Hammer ct al. 
(2009]) for more details about the comparison between 
simulations and observations. This comparison suggests 
that at least 33% of z^0.6 galaxies (50% of galaxies not 
relaxed morpho-kinematically according to the morpho- 
kinematic classifications x the ratio of secured models, 
which is 18/27=66%) are involved in major mergers with 
a median mass ratio ~3. 



4. DERIVING THE MERGER RATE OUT TO Z~1.5 FROM 
MORPHO-KINEMATIC OBSERVATIONS 

4.1. Methodology 

All the merger phases that are expected to 
lead to a significant SFR enhancement, i.e., the 
pre-fusion, fusion/post- fusion, and relaxation phases 
(|Larson fc Tinsle"vlll978|) . are homogeneously populated 
as a function of time by the IMAGES-CDFS galaxies, as 
evidenced by Fig. [4] In this figure, SFRs were homog- 
enized by dividing the observed SFR by the estimated 
mass of gas in each galaxy and multiplying each derived 
number by the average mass of gas in the sample. The 
time sequence reproduces quite well expectations from 
modeling the star formation during a typical merger (see 
insert). This shows that morpho-kinematic data such as 
those provided by the IMAGES survey are sensible to all 
the phases of a gas-rich major merger event. 

Amongst the galaxies populating the relaxation phase, 
six of them are rotating, but they are clearly dis- 
tinct of local rotating disks because of their high av- 
erage st ar formation rate, a nd higher gas veloc ity dis- 
persion (|Puech et al.l l2007t lEpinat et ail 120101 ). The 
latter suggests that these galaxies are possibly ob- 
served during a post major merger relaxation phase 
(|Robertson fc BullockJ |200l iMonreal-Ibero et al.l f20lot 
see also Fig. [3] and Sect. 3.2). Most of the star forma- 
tion activity is indeed concentrated at the disk outskirts, 
which suggests an ongoing inside-out growth of their 
disks (jNeichel et al.l I2008t iHammer et all I2009D . The 
time sequence in Fig. [4] empirically accounts for a cos- 
mological distribution of merger parameters, since it is 
constructed using a representative sample. It can thus be 
used to estimate accurately the evolution of the merger 
fraction and merger rate as follows. 

Theory usually gives the merger rate at the time when 
halos start merging, i.e., where galaxies that inhabit 
these halos are still in pair. The epoch when a given 
merging system was st ill a pair is dated usin g individual 
numerical models from IHammer etldl ([2009D . An homo- 
geneous and self-consistent determination of the merger 
rate as a function of lookback time can be estimated by 
simply counting the fraction of galaxies in each phase 
and dividing this number by the corresponding aver- 
age visibility timescale provided by the modeling of IM- 
AGES galaxies. This method assumes that the visibility 
timescales do not evolve with redshift. Such an assump- 
tion is consistent with semi-analytic models, which pre- 
dict almost no evolution of the merger visibility timescale 
as a func tion of lookback time, at least within uncer- 
tainties (jKitzbichler fc White! 12008). Another possible 
caveat is the obser ved evolution of the galaxy gas frac- 
tion with redshift dPuech et al.l [2 010; Da ddi et all 120101: 
iRodrigues et al.ll2012tlErb et al.ll2006l ). which could po- 
tentially result in evol ving timescales de pending of the 
merging phase probed (jLotz et al.ll2010bl ). H owever, hy- 
drodyn amical simulations of binary mergers ()Lotz et al.l 
20101)) predict that unequal mass ratio mergers (which 
dominate the empirical distribution of merger mass ra- 
tios, see Sect. 3.2) generally show no evolution of the 
pair visibility timescale with gas fraction. 

4.2. The fraction of mergers out to 
4.2.1. Pre-fusion and post-fusion phases 
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spective grade given to the different modeis. It was 
found that the pre-fusion phase corresponds on average 
to z~0.72±g;^ (i.e., 0.7±0.2 Gyr earlier than z=0.6), 
while the post-fusion phase corresponds on average to 



Figure 4. Time sensitivity of morpho-kinematic observations to the major merger process. Galaxies are plotted at the time at which 
individual numerical models give the best-fit to observations. Full symbols represent galaxies whose morpho-kinematics compared best 
with those of their models. Merging galaxies were classified into three different classes: the pre-fusion (blue symbols), post-fusion (red), and 
relaxation phases (green). The pre-fusion phase corresponds to major mergers where both progenitors can be visually identified, while the 
post-fusion phase corresponds to those where the fusion is too advanced to reliably distinguish them. The dash line is not a fit but a simple 
visual guide through the points. The median uncertainty is indicated is the upper-left corne r. The upper-right inset shows typical similar 
curves as reproduced from hydrodynamical simulations of major merger (Cox ct al. 2006). Surrounding stamps show 1MAGES-CDFS 
HST/ACS B-V-z color images with the approximate location of the kinematic axes (dash white lines) and velocity dispersion peaks (white 
dash boxes) superimposed. These stamps are numbered in increasing time sequence: (1) J033224.60-274428.1, (2) J033220. 48-275143. 9, (3) 
J033239. 72-275154.7, (4) J033227.07-274404.7, (5) J033230.43-275304.0, (6) J033245. 11-274724.0, (7) J033237.54-274838.9, (8) J033241.88- 
27 4853.9, and (9) J033 212. 39-274353. 6. Their morpho-kinematic observations are shown in Fig. [6] as well as the result of their modelling 
by I Ham mer ct al. (2009)), see Appendix A. 

Depending on whether or not the two progenitor can- 
didates potentially involved in the merger could be iden- 
tified visually in the HST/ACS images (e.g., when two 
candidates for the progenitor nuclei could be identified), 
the secured major merger candidates were split into pre- 
fusion (8 gal.) or post-fusion phases (10 gal.). Amongst 
the objects that could not be secured as major merger 
candidates using numerical models, three were found to 
correspond to the pre-fusion phase, since they reveal dis- 
tinct potential progenitors. They were therefore added 
to the pre-fusion phase. We attribute the fact that their 
numerical models did not succeed in reproducing their 
morpho-kinematical properties with sufficient precision 
to the limited exploration of the parameter space. 

The pre- and post-fusion phases as observed at z~0.6 
correspond to galaxies that were on average in pairs (and 
before the first passage) at earlier epochs. To date these 
epochs, we simply averaged the time elapsed since the 
beginning of the merger model for each of the 18 secured 
major merger candidates, and converted these numbers 
into a redshift. All numbers were weighted by the re- 



z~1.12±°;^ (i.e., 2.5±0.4 Gyr earlier than z=0.6). Un- 
certainties were estimated using bootstrap resampling. 

The fraction of mergers in emission-line galaxies at 
z=0.4-0.75 for the pre-fusion phase is 0.5x11/33 = 16.7% 
(with the factor 0.5 accounting for the fact that only 
one galaxy in the pair is part of the sample), while the 
fraction of mergers in emission-line galaxies for the post- 
fusion phase is directly given by 10/33 = 30.3%. Since 
60% of interm ediate-mass galaxies at z^0.6 are emission - 
line galaxies (jHammer et all 119971; iMignoli et all 120051) . 
this translates into a total fraction of mergers of 10.0% 
and 18.1% for the pre- and post-fusion phases, respec- 
tively. This assumes that major mergers produce a 
star formation activity that results in EWo([OII]) > 
15A However, it has been shown that major merg- 
ers can lead to m odest SFR enhancement over relatively 
short timescales (|Di Matteo et al.H2008l: ICox et al.|[2008l: 
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iRobaina et~aLll2009l:IJogee et al.ll2009D . This assumption 
can therefore lead to possible systematic effects that are 
considered separately (see Sect. 4.4). Random uncer- 
tainties due to statistical fluctuations in the respective 
sub-samples (see Tab. [1} are derived accurately using a 
beta distribution generator as an est imator of confi dence 
intervals for the binomial statistics (|Cameronll20f ID . 

The f8 secured pre- and post-fusion major merger 
models were used to determine the fusion time for each 
galaxy. The fusion time was defined as the interval of 
time between the beginning of the simulation (i.e., when 
the two progenitors are in pair and before the first pas- 
sage) and the time at which the two nuclei merge accord- 
ing to the best-fit model. These individual timescales 
were averaged to mitigate uncertainties due to statistical 
fluctuations in the respective sub-samples. The averaged 
pre- fusion phase is found to be 1.8±0.1 Gyr long, which 
gives an estimate of the pre-fusion visibility timescale in 
the IMAGES-CDFS sample. The uncertainty was de- 
rived by bootstrap resampling of the grade- weighted dis- 
tribution of individual timescales. We assume the same 
visibility timescale for the post-fusion phase, since both 
appear to be similar in the IMAGES sample, according 
to Fig. m 

We compare this estimate with hy drodynamical sim - 
ulations including radiation transfer ()Lotz et al.l [2008). 
We selected the G2 and G3 models as being the closest 
to the IMAGES-CDFS galaxies in terms of baryonic mass 
and gas fractions. Indeed, the most important parame- 
ters driving the result of a major merger are the baryonic 
mass (as a proxy for the tot al mass inside the opti cal ra- 
dius) and the gas fraction (Hopkin s et al.l 1200 9c) . For 
these two models, the timescale from the beginning of 
the interaction to the fusion of the two nuclei is found to 
be 1.24 and 2.44 Gyr respectively, with an average —1.84 
Gyr, in good agreement with our estimate of 1.8±0.1 
Gyr. We also compared the pre-fusion visibility timescale 
with an average pair visibility timescale for large projec- 
tion distances (i.e., Rnm-i < 100 kpc ./i -1 ) derived from 
similar simulations (jLotz et alJl2010l fbh. which is found 
to be 1.81±0.06 Gyr. Taking into account predictions 
of mass ratios from galaxy formation models slightly re - 
duces the average value to 1.5-1.6 Gyr (|Lotz et al.ll201~Tl ). 
The pre-fusion phase includes galaxies that are too close 
to be selected as pairs by usual criteria (see also Sect. 
4.3). It is therefore expected that the pre-fusion visibil- 
ity timescale in the IMAGES-CDFS sample is slightly 
larger than the average pair visibility timescale. 

4.2.2. Relaxation phase 

The six rotating spiral disks in the IMAGES-CDFS can 
be used to estimate the fraction of mergers at even ear- 
lier epochs (see Sect. 3.2 and 4.1). Their mass doubling 
timescale T s f r is provide d by the ratio betwe en their to- 
tal SFR and stellar mass (jNeichel et alJl2008t ). Note that 
this timescale is independent of the IMF. One of them 
shows a very large T s f r (~7 Gyr) and was therefore con- 
sidered as an interloper and discarded in the following. 
Assuming that these objects are merger remnants, one 
can attribute most of their star formation activity to the 

2 The post-fusion duration was not measured in these simula- 
tions but set to an arbitrary time after the coalescence of the nuclei, 
which prevents any further comparison. 



formation of the disk (see Sect. 4.1 and 6.2). Therefore, 
this timescale can be used to roughly estimate the aver- 
age time elapsed from fusion to the observation epoch of 
these galaxies. On average, these galaxies were therefore 
observed as pairs at a lookback time 1.8+1.8 (pre-fusion 
+ post-fusion phases) + 0.5x(< T sfr >-1.8) = 3.9±0.3 
Gyr, with < T s f r >=2.3±0.3 Gyr (where the error-bar 
was estimated using bootstrap resampling). This corre- 
sponds to z~1.55ig -J}. The average relaxation timescale 
is found to be 5.5-1^8-1.8=1.9 Gyr (see Fig. 0}. 

Amongst the nine non- or semi-relaxed galaxies in the 
IMAGES-CDFS sample that could not be securely iden- 
tified as major mergers, six were found to be in the post- 
merger phase according to their best models. They were 
all modeled quite at the end of the 3 Gyr duration of the 
simulations, which suggests that the morpho-kinematics 
of these galaxies might not have been identified as re- 
sulting from a major merger due to the limited duration 
of the simulations. Amongst them, four show evidence 
for rotation in their velocity fields and were classified as 
perturbed rotators, which is consistent with this inter- 
pretation. We therefore included them as part of the 
relaxation phase (see red points in the relaxation phase 
in Fig. U). This leads to a fraction of mergers at z— 1.6 
of 20.0% (see Tab. Q}. 

4.3. Merger rates 

The merger rate r [Gyr -1 ] corresponding to the three 
merger phases considered above (i.e., the pre-fusion, 
post-fusion, and relaxation phases) can be estimated us- 
ing the ratio between the corresponding merger fraction 
and visibility timescale (see Tab. [T]). Error- bars due to 
random fluctuations in the sub-samples were estimated 
by propagating the random uncertainties on the fraction 
of mergers and on the visibility timescale. The derived 
merger rates and their associated uncertainties are shown 
in Fig. [5j 

For selection criteria similar to those of the IMAGES- 
CDFS sample, a fraction of gal axies in pai r —2- 
6% at z— 0.6 is gene r ally found (IRawat et al.l 120081 : 
IKartaltepe et al.l 120071; IBundv et al.l 120091) . In the 
present study, the fraction of galaxies involved in ma- 
jor mergers is found to be 10±3% at z— 0.7. Such a 
fraction is larger than the fraction of galaxies in pairs 
since the pair technique only picks up progenitors when 
they are approaching for the first time and do not con- 
sider progenitors between the first and second passage 
and/or when they are getting too close one from each 
other, while morpho-kinematic data are sensitive to all 
phases. This is reflected in the visibility timescale as- 
sociated with pair fractions, which is often assumed to 
be —0.5 Gyr from analytic arguments or independent 
simulations (i.e., which do not necessarily refl ect the ex- 
act p roperties of the observations; see also lLotz et all 
1201 1|) . Using such a timescale and a fraction of galax- 
ies in pair —2-6% indeed leads to a merger rate ~4-12% 
Gyr -1 , which is consistent with our z— 0.7 estimate of 
5.5%. When no lower limit to the distance between the 
two progenitors is imposed, the pair technique gives very 
similar results, wit h —9% of galaxie s in pair at z— 0.75, 
and -15% at z-1 ([Ryan et al.ll2008D . 

The morphological technique provides estimates of the 
fraction of morphologically-disturbed mergers of 7-11% 



Galactic disks do not need to survive in A-CDM 



9 



Table 1 

Galaxy merger rate derived from the IMAGES-CDFS survey. From left to right: Merger fraction, merger rate, and redshifts sampled by 
the IMAGES-CDFS survey. Quoted uncertainties correspond to 1-cr random and systematic error-bars respectively. 



Merger phase 


Fraction of mergers 


(%) 


Merger rate (%.Gyr L ) 




Redshift 


Pre-fusion 


i n n^^ 1 ^ +3.3 
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c cr + 3.54 
°-°-1.8- 


-1.9 
-2.0 


2 


- 72+ u - lU 


Post-fusion 


IS 1+6-7+3.3 
t»-t_3. 8 _5.5 




10.1+" 


+1.8 
-3.3 


Z 


_ , , 9 +0.13 
— LJ^-O.IO 


Relaxation 


on n+6-8+3-3 
—4.0—9.1 






+1.8 
-5.1 


Z 


_ i „+0.15 



at z~0.7, 11-17% at z~l.l, and 27±11% at z-1.4 for 
similar mass ranges an d depending and the method used 
(CAS /GM system, see iConselice et al.1 120081: lLotz et al.1 
2001 [Conselice et all 120091 or vis ual inspection, see 
Bridge et al.H20ia IJogee et al.H2009[ ). Within uncertain- 
ties, this is consistent with the fraction of mergers found 
using morpho- kinematic data (see Tab. [TJ. 

We conclude that alternative techniques generally give 
consistent esti mates, although wit h a resulting large 
range of values ( Hopki ns et al.l l2009). which reflects sys- 
tematic effects due to their limited sensitivity to the 
merging process and di fferent sample selection effects 
(see also lLotz et al.|[201lh . 

4.4. Systematic uncertainties 

A first possible cause of systematic effects is linked to 
the selection criteria of the IMAGES-CDFS sample, and 
in particular to the criterion E Wq ([0/7]) > 15 A. Indeed, 
we assumed that all quiescent (EWq([OII]) < 15 A) 
galaxies are not major mergers, whereas a significant 
fraction of them (25%) were found to h ave peculiar mor- 
phologies (Dclgado - Serrano et al.|[2010() . which could be 
linked to less active mergers. Indeed, numerical sim- 
ulations have revealed that major mergers do not al- 
ways lead to strong enhancemen ts of the star forma- 
tion rate (|Di Matteo et al.1 12008H , and therefore these 
peculiar quiescent galaxies might also be undergoing 
major merger events. This could potentially rise the 
fraction of mergers by 10.0% since quiescent galaxies 
represent 40% of intermedia te-mass galaxies at z~0.6 
(jDelgado-Serrano et al.1 l2010f ). Since we do not have 
kinematic data for these objects, this could potentially 
affect any of the three merger phase considered in this 
study. We therefore considered that, on average, the 
merger fraction in each phase could be affected by a 
~3.3% systematic effect upward (see Tab. Q]). 

The number of mergers in the pre-fusion phase 
could also be overestimated due to fly-by interlop- 
ers. We estimated that four out of the eleven ob- 
jects in this phase could be possible fly-byes and escape 
a subsequent merger (J033217.62-274257.4, J033219.32- 
274514.0, J033220.48-275143.9, J033238.60-274631.4). 
These objects correspond to cases where the two po- 
tential progenitors are well separated and photometric 
redshifts are consistent, but no reliable spectroscopic in- 
formation was available to confirm whether or not the 
two galaxies were indeed gravitationally bound. This 
can potentially have an impact on the fraction of merg- 
ers in the pre-fusion phase by 3.7% downward, and the 
corresponding merger rate by 2.0% (see Tab. [T]). 

Six galaxies in the IMAGES-CDFS are consisten t with 
expectations from clumpy galaxies (Pucch 2010), i.e., 
gas-rich Toomre-unstable disks. These six galaxies are 
found to be in the post-fusion phase (3 out of 6). Since 



the numerical merger models do not consider small-scale 
perturbations such as clumps, these galaxies might cor- 
respond to Toomre-unstable disks driven by internal in- 
stabilities, rather than major mergers. The post-fusion 
fraction of merger and merger rate might therefore be af- 
fected by a 5.5% and 3.3% effect downward, respectively 
(see Tab. [TJ. The impact on the pre-fusion and relax- 
ation phases is marginal (one object) compared to other 
uncertainties, so we did not consider it further. 

Regarding the relaxation phase, we added six objects 
corresponding to the post-fusion phase but which were 
not safely identified as mergers. This could be due to in- 
sufficiently long evolution times, so that the objects could 
not be safely identified, or to a too limited exploration of 
the parameter space. Considering these six galaxies as 
possible interlopers for the relaxation phase, this would 
impact the corresponding fraction of merger and merger 
rate by a 9.1% and 5.1% effect downward, respectively 
(see Tab. [Q). 

Systematic uncertainties are shown as a grey region in 
Fig. [5] 

5. COMPARISON WITH SEMI-EMPIRICAL MODELS 

Instead of trying to predict galaxy properties all the 
way from initial conditions given by the CMB, semi- 
empirical models by-pass the unknown coupling between 
dark and baryonic matter using empirical relations, such 
as gas fraction or stellar mass functions as a function 
of redshifts. Such an approach allows us to overcome 
several issues encountered with semi-analytic models, 
whic h can lead to significan tly underestimate the merger 
rate (Hop kins et all l2010bl). Th e semi-empirical model 
used here (|Hopkins et al.l [2009) was shown to repro- 
duce a number a local galaxy properties, such as their 
bulge-to-disk flux ratio distribution, Tully-Fisher rela- 
tion, or luminosity fun ction of luminous infrared galaxies 
([Hopkins et al.ll2010d ). 

5.1. Model parameters 

We used the iHopkins et al.l (|2009ft semi-empirical 
mode0 to predict the merger rate in a given range 
of baryonic mass [M™ m (z),M™ ax (z)}, gas fraction 
[/go" l ( z ): and baryonic mass ratio between 

the two progenitors \/j,™ m (z),iJ,™ ax ( z)] as a function of 
redshift z. Indeed, iHopkins et al.l (|2009l lcl) emphasized 
that the most appropriate mass would be the tightly 
bound material that survives stripping to strongly per- 
turb the primary, which is found to be the baryonic 
plus tightly bound dark matter mass within a few 
disk scale lengths. The best obs ervable proxy f or thi s 
quantity is the baryonic mass ([Hopkins et al.l [2009), 

3 available at https://www.cfa.harvard.edu/~phopkins/Site/ 
mergercalc.html as an IDL code. 
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Figure 5. Galaxy merger rate out to z~2. The black points with horizontal and vertical error bars show the observational estimates of the 
merger rate r from the IMAGES-CDFS sample. The grey-shaded region shows the impact of possible systema tic effects of the obser vational 
estimation of the merger rate (see Tab. [TJ. The red curve shows the prediction of the semi-empirical model ( Hopkins ct al. 2009), for the 
range of baryonic mass and mass fraction observed. The shaded red region shows the random uncertainty associated to the theoretical 
prediction (see Sect. 5.2). The dashed-red line shows the impact of systematic uncertainties in stellar mass stellar population models on 
model predictions (see Sect. 5.2). The upper-left inset shows the cumulative merger rate R as a function of lookback time (i.e., from the 
full-red curve of the main panel). 



which was therefore used in this study. In most stud- 
ies, the galaxy merger rate is estimated in a fixed 
range of stellar or baryonic mass, constant with red- 
shift, i.e., [M b mm (z),M b max (z)} = [M b min ,M b max }. Here, 
we need to estimate the proper range of baryonic mass 
[MJ"™(z),M 6 mra (z)] at all z in order to sample the pro- 
genitors and descendants of the IMAGES-CDFS galax- 
ies. 



Given [M b min (z),M b max {z)} 1 we determined the corre- 
sponding range of halo mass [M™ m (z),M™ ax (z)] us- 
ing the relationship between baryonic and halo m asses 
M b - M h as provided by the iHopkins et~ail (|2009[ ) halo 
occupation model (H O D) us ing the mass function from 
iPerez-Gonzalez et al.l (|2009f) . Using alternative mass 
functions change the merger rate by less than 1%. We 
then estimate the ave rage mass of the progenito r ha- 
los at z + dz following iNeistein et al.1 (j2006l |2008| ) . i.e., 
[M™ m (z + dz),M™ ax (z + dz)}, which is then converted 
back to [M b mm (z + dz),M b max {z + dz)} using the M b -M h 
relationship. The HOD model is used to estimate the 
galaxy merger rate at z + dz in the updated range of 
baryonic mass. This ensures that the derived galaxy 
merger rate is integrated over the population of the 
progenitors of all galaxies within [M™ m (0),M™ ax (0)] 
at z = only, and not over all the mass spectrum, 
which would overestimate the merger rate compared 
to estimates from observations. In practice, we set 
[M™ m (0)/M Q , Af 6 maa: (0)/M Q ]~[10 10 - 4 ,10 11 - 1 ] in order 
to sample the range [M™"(z = 0.6)/M Q ,M™ aa: (z = 



O.6)/M ] = [lO lo - 2 ,lO n -°], which corresponds to the bary- 
onic mass range of the IMAGES-CDFS sample (see Sect. 
2). We tested that following the mean progenitor mass 
instead of the average progenitor mass between the up- 
per and lower limit of the mass range does not change 
the resulting merger rate significantly. 

To test the consistency of the adopted HOD model 
with the IMAGES sample, we checked that, ac- 
cording to the HOD model, the range of baryonic 
mass sampled at z=0.6 corresponds to a range in 
log(M steHar /M©)=[10 10 - ,10 11 - ], which also matches 
the range of observed stellar mass in the IMAGES sample 
(sec Sect. 2.1). At z=0, the HOD model predicts that the 
z=0.6 progenitors as defined by their range of baryonic 
mass as above, will end up in a stellar mass range of 
log(M s t e/ter /M )=[lO lo - 3 ,lO u - 1 ] ) corresponding to the 
baryonic mass range of log (Mb/M )=[lO lo - 4 ,lO 11 - 1 ] as 
quoted above. The evolution in stellar mass as predicted 
by the HOD model is therefore ~0. 1-0.3 dex. 

The median gas fraction in the progenitors of 
the IMAGES-CDFS s ample were estimate d using 
their observed SFRs (jHammer et al.l 120091 ). We 
used this to parameterize the expected gas frac- 
tion in the IMAGES-CDFS progenitors and de- 
scendants, using / sas ={15%,32%,50%,75%,90%} at 
^={0,0.6,0.72,1.12,1.55}. Estimates at z=0 and 0.6 were 
derived using local HI measurements and in version of the 
Schm i dt-Kennicutt relation, respectively (jPuech et al.l 
120101: IRodrigues et al.l I2012D . Estimates at z=0.72, 
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1.12, and 1.55 were derived from the estimated median 
gas fractions in the pre-fu sion, post-fusion, and relax- 
ation phase, respectively (jHammer "eFaD [2009h . The 
HOD model predicts that the progenitors of z— 0.6 on- 
going mergers in the IMAGES-CDFS sample are galax- 
ies with stellar mass in the range log (M s t e iiar/M@) ~ 
[10 9 2 ,10 10 - 4 ]. Gas fractions of the order of -90% in 
such a range of mass at z— 1.5 is consistent with observa- 
tional estimates of the gas fraction in more z—1.5 mas- 
sive galaxies (see discussion in Sect. 6.2). We added 
a ±l-cr scatter of 0.15 dex, which app ears to repre- 
sent we ll the observed scatte r at all z (H opkins et al.l 
l2009H bl: iRodrigues et al.l 12012ft . At given z, the range 
of gas fraction [f^ l J l (z),f^ x (z)] over which the merger 
rate was integrated is interpolated from the f gas (z) re- 
lationship as defined above, with an added ±2-cr scat- 
ter. We explored the impact of a larger ±3-cr scatter 
with no qualitative change in the results. Finally, we 
set \uF in (z),u™ ax (z)] = \ur n ,u r ax ] = [0 25, 1.0] as fo und 
in the IMAGES-CDFS sample (jHammer et al.|[2009h . 

5.2. Comparing observations and theory 

Figure [5] compares the derived estimates of the merger 
rate with state-of-the -art semi-empirical A-CDM models 
(Hopk ins et alj 120091 ) . This theoretical model is found 
to be in remarkable agreement with the observed esti- 
mates without any particular fine-tuning. At the red- 
shifts and masses probed by the IMAGES-CDFS sam- 
ple, semi-empirical (HOD-based) models can predict th e 
merger rate within a factor — 1.5 (Hopkins et al. 2010b). 
HOD models are in relatively good agreement one with 
each other because they rely on observational constraints 
(e.g., the galaxy stellar mass function) that are well 
constrained at z<l a nd for stellar masses around M* 
(Hop kins et al.l l2010bf l . As stated, above, changing the 
galaxy stellar mass function changes the resulting merger 
rate by less than — 1%. The dominant source of ran- 
dom uncertainty is therefore due to variations from one 
(HOD) model to another, which is illustrated in Fig. [5] 
as a red region. This is the same order of magnitude that 
random uncertainties associated to the observed merger 
rate (see Tab. [I}. This implies that random uncertainties 
limit the comparison between observations and theory at 
a level of a factor —2. 

In this comparison, we adopted a diet Salpeter IMF 
for all observables at any epoch. Several studies sug- 
gest that the IMF might b ecome top-heavy un der ex- 
treme starburst conditions (We idner et al.l f2010h . This 
could result in an evolution of the IMF with redshift, 
which would impact the estimation of stellar mass by 
up to — -0.18 dex (rising almost linearly with red- 
shift from z=0 to z— 1.5) compared t o estimates at z=0 
(jvan Dokkum fc van der Marell 120071 ). This could result 
in a systematic shift of the merger rate with redshift. 
Several systematic effects associated with the stellar pop- 
ulation models can also affect the determination of the 
stellar mass. The most significant effect is the influence 
of the TP-AGB phase, which omission in stellar popula- 
tion models as those used in the present st udy can lead to 
overe stimate the stellar mass by 0.14 dex (|Pozzetti et al.1 
120071) . In addition, using simplified prescriptions for 
deriving the stel lar-mass instead of full SED modelling 
(jBell et al.ll2003[ ). when used at high redshifts, can also 



potentially lead to overestimate stellar m asses by 0.2 dex 
(|Puech et al.1l2008t iG^hazzi fc Bellll2009h . 

In order to bracket the upper limit of the predicted 
merger rate at all redshifts, we combined quadratically 
the systemic uncertainty associated with a possible evo- 
lution of the IMF, and the effects of stellar population 
models, since they are independent. This results in a 
total systematic uncertainty of -0.27 dex on the bary- 
onic mass. In Fig. [31 we plotted theoretical predictions 
with [Ml nm (z),M™ ax (z)]-0.27, using a dashed-red line. 
In doing so, [f^T( z )Jgas X ( z )] was updated accordingly 
(upward). All these systematic uncertainties associated 
to the estimation of the stellar mass tend to decrease 
the galaxy merger rate at all redshifts. However, the 
theoretical predictions remain consistent with the obser- 
vational estimates within the corresponding range of ob- 
servational systematic uncertainties. Quantitatively, this 
could decrease the integrated merger rate by —20%. 

We conclude that the major merger rate is deter- 
mined from both theory and observations within a fac- 
tor 2-3, accounting for random and systematic uncer- 
tainties. Such an agreement between the merger rate 
determined using morpho-kinematic data and theoreti- 
cal predictions without any particular fine-tuning, is re- 
markable. This is an improvement of a factor —3 com- 
pared t o previous compariso ns between observations and 
theory ([Hopkins et al.ll2009f ). 

6. DISCUSSION 

6.1. One third of z— 0.6 galaxies are on-going major 
mergers 

Using morpho-kinematic data of z— 0.6 galaxies from 
the IMAGES-CDFS survey, we derived the evolution of 
the major merger rate out to z— 1.5, i.e., nine billion 
years ago. Using comparisons of the spatially-resolved 
morpho-kinematic properties between observations and 
a grid of hydrodynamical simulations of galaxy mergers, 
we constrain their time evolution, which allowed us to 
accurately date back when they were still in pair. In this 
comparison, the estimate of the merger rate at z— 1.5 
is much more uncertain since it relies on the relaxation 
phase during which it is more difficult and / or degenerate 
to interpret morpho-kinematic features as merger relics. 
However, it makes sense to interpret such galaxies as 
post-mergers in the relaxation phase as detailed above, 
and Fig. [3] & 0] indeed suggest a plausible causal link 
between distant merging galaxies and such dynamically 
hot rotating disks. We found a remarkable agreement 
between these new observational estimates and predic- 
tions from state-of-the-art semi-empirical A-CDM mod- 
els within a factor 2-3, accounting for both random and 
systematic uncertainties. 

As a test of consistency, the A-CDM model predicts a 
large fraction of peculiar galaxies at z=[0.4-0.75], which 
is found to be < r > (i.e., the average merger rate at 
z=0.65, -0.06 Gyr" 1 , see Fig. [5j) X r (i.e., the total 
duration of the average major merger at z=0.65, —5.5 
Gyr, see Fig. [4| — 33%. This corresponds exactly to 
the fraction of observed galaxies in which the peculiar 
morphologies and/or kinematics can be reproduced by a 
major merger model (see Sect. 3.2). The A-CDM model 
therefore predicts that —33% of the progenitors of local 
disks were destroyed over the past 9 Gyr during major 
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(i.e., baryonic mass ratio >1:4) merger events, as ob- 
served. 

The occurrence of 33% of gas-rich major mergers in 
z~0.6 galaxies is also consistent with an increase of 
the scatter of the Tully-Fisher relation at high redshifts 
(|Flores et alJl200a iKassin et al.ll2007t iPuech et all[200l 
iCovington etalll2010t iPuech et al.l[2010h 0, and the ob- 
served random walk in the specific angular momentum of 
z~0.6 galaxies compared to local galaxies (|Puech et al.1 
l2007h . 

6.2. Disk rebuilding of local bulge-dominated spirals 

The agreement found between observations and the- 
oretical predictions (see Fig. [5J evidences that 50% 
(65%) of local galaxies with stellar mass in the range 
1O 1O ' 3 -1O 111 M0 (see Sect. 5.1) have undergone a ma- 
jor merger (defined here a s baryonic mas s ratio >0.25) 
since z=1.5 (z=2.0). The iHopkins et all (|2009f ) model 
takes into account the specific physics involved in gas- 
rich major mergers, which was shown to reproduce the 
observed loc al relation between b ulge-to-disk ratio and 
stellar mass (|Hopkins et al.l 2009. b). This suggests that 
the time at which local galaxies have undergone their last 
major merger could be an important factor in determin- 
ing how they acquired their internal structure (B/T ratio 
and angular momentum). 

Depending on the gas fraction at fusion time, ma- 
jor mergers are well-known to result in ellipticals 
(through dry mergers), but also in spirals (through 
gas-dominated mergers). Indeed, hydrodynamical sim- 
ulations of isolated merging disk galaxies have shown 
that with a fraction of gas at the fusion time larger 
than 50%, a significant disk can reform in the rem- 
nant (iSpringel & HernauistJl2005HRobertson et al.ll2006t 
Hop kins et all l2009d) . Note that the required fraction 
could be smaller with sig nificant cosmological gas accre- 
tion ([Brooks et all 120091 ) , but we will adopt 50% as a 
conservative limit. 

The median gas fraction at z=0.6 is fou nd to be ^30% 
(jPuech et al.ll2010HRodrigues et al.ll2012fl . It is therefore 
likely that major mergers occurring between z~0.6 and 
z=0 (which involved 20% of local galaxies according to 
Fig. [5]) were relatively dry, probably resulting in galaxies 
with significant bulges. 

The major mergers sampled by the IMAGES-CDFS 
sample are those occurring between z^l.5 and z~0.6 
(see Tab. []}. Fig. [5] shows that 50% of local galaxies 
were involved in a major merger since z =1.5, of which 
30% occurred between z=0.6 and z=1.5. Hammer et al. 

4 It has been recently claimed that the distant stellar mass TF 
relation might have a scatter that is actually similar to the local 
relation (Miller et al. 2011). The reduced scatter has been at- 
tributed to the long integration times that would better sample 
the plateau of the rotation curve combined to an optimal velocity 
extraction method, providing a more accurate circular velocity es- 
timation. We note that exposure times in the IMAGES survey are 
comparable or even larger that those used in this study. Moreover, 
IFU data provide additional element of resolutions in the direc- 
tion perpendicular to the galaxy major axis, which provide more 
information to accurately deconvolve the line-of-sight velocity mea- 
surement and estimate the rotation velocity. This is essential when 
dealing with galaxies with complex VFs, since the kinematic and 
morphological axes can be significantly misaligned, or when decon- 
volving the measurement from seeing effects. It is beyond the scope 
of this paper to provide a definitive assessment on the evolution of 
the scatter of the TF relation with redshift. 



(2009) have extrapolated the gas fraction in the progen- 
itors of these mergers, and found gas fractions well in 
excess of the required value of 50% to reform a disk 
(~90%, see Sect. 5.1). Several me asurements of the 
gas fraction were att empted at z~1.5 (jDaddi et alll2010t 
iWeinzirl et all 120111 ). Most measurements involve rela- 
tively massive galaxies (~M ste ;/ ar = 10 ia7 — 10 111 M Q ) 
and are limited to their half light radius. However, the 
HOD model predicts that the progenitors of the on-going 
z~ 0.6 major mergers at z=1.5 sh ould be in a lowe r range 
of stellar mass -10 9 ' 2 -10 10 - 4 M Q . iRodrigues et all (|2012f ) 
homogenized gas estimates from the literature to param- 
eterize the evolution with redshift and stellar mass of 
the gas fraction within the full gas radius. Extrapolating 
their gas fraction-stellar mass relation at z^l.5 in the 
range of mass 1O 9 2 -1O 1O 4 M0, one finds that gas frac- 
tions at a level of ~ 80-100% ar e probably close to av - 
erage conditions (see Fig. 7 of IRodrigues et ail [20121 
The progenitors of all the z^0.6 on-going major mergers 
at z=[0.75-1.5] should therefore be rich enough in gas so 
that a significant disk reform after fusion, as also evi- 
denced b y observations of young dust-enshrouded disks 
at z~0.6 (jHammer et al.ll2009bUPuech et aLll2009D . 

We emphasize that we have considered here only the 
"average" results of major mergers since z~2 and such 
statements have to be taken at the lowest order only, 
since one expects significant scatter in the B /T of a ma- 
jor merger remnant as a fun ction of orbital parameters , 
gas fraction, and feed back ([Hopkins et all 120081 120091: 
iGovernato et al.ir2009D . We also emphasize that we con- 
sidered here only a specific range of mass and that the 
impact of more modest mergers, i.e., with baryonic mass 
ratio smaller than 0.25 was not considered. 

6.3. Is there a disk survival issue? 

Results of semi-analytic models suggest that the likeli- 
est channel for the formation of disk-dominated galaxies 
is a relatively qu iet merger history ( e.g., see discussion 
and references in IWeinzirl et al.ll2009f ) . Of particular in- 
terest are local "bulgeless galaxies " (which we define 
here as galaxies with B/T<0.2, as generally considered), 
which are often thought to have undergone their last ma- 
jor merger above z~2. The relatively large fraction of 
such galaxies is often considered to be in tension with the 
fraction of major me rgers as inferred from predictions of 
A-CDM models (e.g- lStewart et al. 2008). It was argued 
in Sect. 6.2 that most descendants of IMAGES-CDFS 
galaxies should lie within early- type spirals. We now ex- 
amine whether the fraction of mergers found in Fig. [5] 
is in tension with the observed fraction of local bulgeless 
g alaxies. 

IWeinzirl et all (|2009l ) used Bulge+Bar+Disk decompo- 
sitions on H-band images and found that two-third of lo- 
cal spirals with M ste M ar =lO lo - 1 -lO 11 - 6 M are bulgeless. 
Their determination maximizes the fraction of bulgeless 
galaxies since they removed the contribution of the bar 
to the central light and followed the common assumption 
that bars are pure products of secular evolution and not 
merging (see discussion in I Athanassoulal 120081 ). Since 
the fraction of spirals amongst local galaxies is ~70% 
(see App. B), this would translate into a maximal frac- 
tio n of bulgeless g alaxie s in the local Universe ^46%. 
As IWeinzirl et all (|2009l ) argued, their sample appears 
to be representative of the B-band luminosity function 
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at M s teiiar > 10 10 ' 1 M (see their Fig. 2). However, their 
Fig. 3 evidences that their sample can be representative 
of stellar mass function only above M ste ;/ ar >10 10 6 Mq. 
In this range of mass, the maximal fraction of local bul- 
geless is found to be 61%, or 43% of local galaxies. As- 
suming that such galaxies cannot be formed by mergers, 
this would therefore limit the fraction of major merger 
in local galaxies since z=2 to 57%. We re-derived Fig. 
[5] for the range of mass M steH ar=10 10 - 6 -10 1:L6 M Q . In- 
deed, the HOD model used in the present study predicts 
that on-going major merger in the z^0.6 IMAGES-CDFS 
sample should end up in a stellar mass range of 10 10 ' 3 - 
1O 1:L1 M0, which is shifted to smaller m ass relative to 
the range 

10 io.6_ 10 ii.6 Mq in which the |Weinzirl et ail 
(2009) sample is representative. In this higher range of 
mass, we found that the fraction of local galaxies that 
have undergone a major merger since z=2 decreases from 
65% to 53%. Within uncertainties, this matches quite 
well the maximal fra ction of bulgeless galaxies found by 
iWeinzirl et all (|2009t) in the same range of mass (57%). 

We conclude that there is no tension between the frac- 
tion of mergers inferred from morpho-kinematic data and 
the fraction of local bulgeless galaxies. In other words, 
there is no need for disks to survive major mergers in the 
A-CDM paradigm, at least at z<2. 

6.4. Disk rebuilding in the cosmological context 

The results of this paper point toward a re-building 
of half of local disks during th e past 9.4 Gyr aft er 
their last gas-rich major merger iHovkins et 
\Hammer et al\\200A : \Stewart et al\\200<\ ). This does not 
imply that the evolution of the star formation and stel- 
lar mass densities are primarily driven by such processes 
over th e same period (jRobaina et al.ll2009l : lHopkins et al.l 
l2010d ). These results are consistent with recent de- 
velopments of the A-CDM model at earlier lookback 
times (i.e., z>1.5), according to which galaxy formation 
would be mainly driven by streams of cold gas, rather 
than by major mergers (iForster-Schreiber et al.l 120091 : 
lAgertz et al.fl2009l : IWeinzirl et al.ll2011l ). Indeed, cosmo- 
logical simulations predict that cold flows a re strongly 
suppressed below z^l-1.5 (jKeres et al.1 12009). Such an 
epoch might correspond to a transition between two dif- 
ferent mechanisms driving galaxy formation, i.e., cold 
gas streams and gas-rich major mergers, as suggested 
by the strong evolution of the fraction and properties of 
clumpy galaxies as a function of redshift (Puechl 120101 : 
IContini et al.l l2012b[ ). In such simulations, halos with 
a relatively quiet merger history are usually selected as 
the best probes fo r the formation of local thin disks (e.g., 
iBrook et al.ll2012[ ). The present results suggest that gas- 
rich major mergers could be be an equally plausible way 
of forming large thin disks. Further work will be needed 
to quantify what is the main reservoir of gas from which 
such thin disks can reform. Potential sources of gas 
include gas expelled du ring the merg er and which falls 
back onto the remnant flBarnesI 120021), gas coming from 
external gas accretion (jBrooks et al.l I2009T) . and/or gas 
recyc led through stellar feedback (jMartig &: Bournaudl 
120101) . This emphasizes the need to further explore how 
such disks can (re-)for m in halos with a mo re recently 
active merger history (|Governato et al.l 12009). Interest- 
ingly, recent cosmological simulations conducted with 
state-of-the-art hydrodynamical codes reveal cold gas 



disks that have larger extent and specific angular mo- 
mentum compared to d isks formed in former SPH cos- 
mologicial simulations (|Torrev et al.l 120111 : iKeres et al.l 
120111 ). These simulations show that massive halos (i.e., 
Mhaio > 10 11 Mq) have significantly l arger SFRs at z=0 - 
2, compared to previous simulations ()Keres "elalll20ll . 

7. CONCLUSION 

Using morpho-kinematic data of z=[0.4-0.75] galaxies 
from the IMAGES-CDFS survey, we derived the evolu- 
tion of the major merger rate out to z^l.5, i.e., nine 
billion years ago. Using comparisons of the spatially- 
resolved morpho-kinematic properties between observa- 
tions and a grid of hydrodynamical simulations of galaxy 
mergers, we constrain their time evolution, which al- 
lowed us to accurately date back when they were still 
in pair. We then compared constrains from observations 
with predictions from semi-empirical models. Our main 
conclusions are as follows. 

1. Morpho-kinematic observations are found to be 
sensitive to all the phase of the merger process, 
from the pre-fusion phase in which galaxy nuclei 
are still separable by visual inspection, to the post- 
fusion and relaxation phases, in which the two pro- 
genitors merge and progressively reach equilibrium. 
Each of the three phases are found to last ~1.8 Gyr. 

2. The major merger rate (defined as baryonic mass 
ratio >0.25) inferred from morpho-kinematic ob- 
servations is found to be evolving from 5. 5%. Gyr - 
to 11.1%. Gyr" 1 between z=0.72 to z=1.55 in the 
range of stellar mass 10 10 -10 11 Mq (diet Salpeter 
IMF), or lO lo - 2 -lO n M in baryonic mass. 

3. A remarkable agreement is found between ob- 
servational es t imate s and predictions from the 
Hopk ins et al.l (|2009f) semi-empirical model, which 
accounts for the specificity of gas-rich mergers at 
high redshifts. The agreement is found to be at 
a level of a factor 2-3 without any particular fine- 
tuning of the model, accounting for both random 
and systematic uncertainties on the observational 
and theoretical sides, which significantly improves 
the accuracy of previous comparisons. 

4. The A-CDM model predicts, in agreement with 
observations, that 33% of z~0.6 intermediate- 
mass galaxies are on-going major mergers. Ob- 
servational signatures of these mergers encompass 
morpho-kinematic disturbances in a large fraction 
of these galaxies at this redshift, and an increase of 
the scatter of scaling relations between mass, ve- 
locity and radius (or angular momentum). 

5. The A-CDM model predicts, in agreement with ob- 
servations, that 20% / 50% / 65% of local galaxies 
in the range of stellar mass 1O 1O 3 -1O 111 M0 were 
involved in a major merger between z=0.6 / 1.5 / 
z=2.0 and z=0. 

6. Most of the major mergers occurring between 
z=0.6 and z=0 were probably relatively dry, result- 
ing in galaxies with a prominent bulge, while most 
of the major mergers occurring between z=1.5 and 
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z=0.6 probably involved large gas fractions, which 
should lead to the rebuilding of a significant disk. 

Theoretical models often suggest that "bulgeless 
galaxies" (B/T<0.2) result from relatively quiet 
merger histories. We found no tension between 
the fraction of local bulgeless galaxies and the frac- 
tion of local galaxies that have undergone a major 
merger since z=2, i.e., there seems to be no "disk 
survival issue" but rather a "disk rebuilding" era, 
at least at z<1.5-2, when the influence of cosmo- 



logical gas accretion predicted by theory and sim- 
ulations would start declining in relatively massive 
haloes. 
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APPENDIX 

EXAMPLES OF MORPHO-KINEMATICS DATA AND ASSOCIATED MAJOR MERGER MODELS IN THE IMAGES-CDFS 

SAMPLE 

We present here the i-band HST/ACS images and kinematic maps for a sub-sample of the IMAGES-CDFS, see Fig. 
O Galaxies are the same that those used to illustrate the SFR vs. time sequence of Fig. [H in the same (i.e., temporal) 
order, from top to bottom. Galaxies with no major merger model are RDs. 

EVOLUTION OF GALAXY MORPHOLOGY BETWEEN Z~0.6 AND Z=0 

We discuss here the evolution of the fraction of morphological types between z~0.6 and z=0 in the range of stellar 
mass which encompasses the descendant of the IMAGES-CDFS galaxies. 

We estimated the fraction of each morphological type in local gala xies as follows. We first considered Morphologically 
Determined Luminosity Functions (MDLF) of local galaxies from iNakamura et al.l (|2003l ). They indeed constructed 
LFs for E/SO, SO/a-Sb, Sbc-Sd, and Im galaxies using ^1500 galaxies from the SDSS usin g the r band lum inosity. 
We first converted their observed r band luminosity into stellar mass using prescriptions from I Bell et al.l (|2003f ) . which 
gives a mean stellar mass/r-band luminosity ratio in the local volume of 3.05 within ±30% of systematic uncertainty 
(diet Salpeter IMF). We then integrated each of the MDLF (once converted into the proper IMF) from 10 10 3 to 
10 111 M©, within which the descendants of the on-going z^0.6 major mergers should end up, as predicted by the HOD 
model (see Sect. 5.1). This results in the fractions of morphological types as listed in Tab. [3] Deriving the fraction 
of morphological type that way has the advantage of avoiding systematics associated with th e representativity of the 
sample, and is limited only by the statistical scatter associated with the visual classification of Nakamu ra et al.l (|2003l ) 
and the sy stematic uncertainty associated to the mean stellar mass-to-light ratio in the local volume as determined by 
IBell et all (pOOl . 

As a consistency check, we also derived the fraction of morphological types i n the local Universe using the catalogue 
of ^14000 visually classified galaxies in the SDSS from Nai r fc Abraham! (|2010f ). We selected in their catalogue galaxies 
in the range of mass 10 la3 -10 11 ' 1 M Q , and with redshift z < 0.0 5. According to their Fig . 1, this select a volume-limited 
sample in the stellar mass range of interest for our purposes. iNair fc Abraha m (2010) emphasized that such a cut in 
z will remove some of the brightest, most massive galaxies in the local volume, but these galaxies are all above the 
range of mass of interest, so this will not impact the morphological fractions in this range. Derived fractions are also 
listed in Tab. [2j It is interesting to compare them with those derived by integrating the MDL Fs in the local volume, 
since b oth derivations rely on completely different sel ection effects. We also listed results from iDelgado-Serrano et all 
(2010), which selected a sample of the catalogue from iFukugita et all (|2007l ). who visually classified a sample of more 
than 2000 galaxies from the SDSS. Their sub-sample was selected to be volume-limited with M ste ;/ ar >1O 1O - 3 M to 
specifically match the descendants of the galaxies sampled by the IMAGES survey. 

The comparison between these three different determinations of the fraction of morphological types in the local 
Universe in the same mass range reveals interesting features. First, the fraction of spirals is found to be in remarkable 
agreement, around ^70%. There are some variations within the fractions of early-type spirals, probably reflecting 
the uncertainty associated to eyeb all classification betwee n SO a nd SOa galaxies. The second interesting feature is the 
smaller fraction of E/S0 found by IDelgado-Serrano et all (|2010l) in conjunction with their hig her fraction of Peculiar 
galaxies compared to other catalogues. This was attributed bv [Delgado- Se rrano et all ((2010) to a number a factors 
that include (1) their systematic use of color maps to identify unusual color gradients, (2) the quantification of the 
misalignment between bulge and disk centroids, and (3) the compactness of galaxies, which were moved to the peculiar 
cl ass if too compact to be decom posed into bulge and disk. Galaxies harboring such features were classified as Peculiar 
bv IDelgado-Serrano et al.1 ((2010), but are probably in the E/S0 class in o ther catalogues that do not take into account 
such a n information. This suggests that ~10% of the E/S0 classified by INakamura et ahl (|2003l ) or INair fc Abraham! 
(2010) arc actually perturbed and not fully relaxed. Such an interpretation is consistent with the fraction (~30%) 
of early-type gal axies showing negative color gradients ("blue-cored" early-types) in the local Universe in the SDSS 
(|Suh et a.l.ll2010h . 

At face value, the fracti on of late- type spirals in Tab. [H i.e., 20-23% might appear in conflict with the fraction of 
bulgeless spirals derived bv lWeinzirl et al.1 (|2009f ). which is found to be ^60% (see Sect. 6.3). We recall that "bulgeless" 
means her e B/T<0.2 and that the correlation between B/T and morphological type has very large scatter (see, e.g., 
Fig. 14 of iWeinzirl et al.l 1200 9) . Selecting galaxies with morphological types in the range Sbc-Sd in the sample of 
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Figure 6. Morpho-kinematic observations for an IMAGES-CDFS subsample and results of their modelling. See text for details. For each 
galaxy, from left to right: HST/ACS i-band image with the GIRAFFE IFU superimposed, GIRAFFE velocity field, and GIRAFFE velocity 
dispersion map (top panel), modelled stel lar density, modelled gas density, simulated velocity field, and simulated velocity dispersion map 
(bottom panel). Sec Hammer ct al. (2009) for details. 
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Figure 7. Continued. 



IWeinzirl et al.l (|2009D leads to a fraction of late-type galaxies ~33% (above M ste ;; ar =lO 1O 6 M0). Since spirals account 
for ^70% of local galaxies (Tab. [5]), this translates into a fraction of ~23% of late- types amongst local galaxies, which 
is fully consistent with fractions listed in Tab. [5] 

We also listed in Tab. [2] the fraction of morphological types found by iDelgado-Serrano et al.l (|2010f ) at z^0.6. Note 
that the fractions of Peculiar galaxies, i. e., those that do not fi t into the Hubble sequence, a re at both redshifts in 
very good agreement with those found bv lvan den Berghl (|2002t l. [Delgado-S errano et al.l (j201Q[ ) used exactly the same 
classification method and strictly equivalent data to classify local and distant galaxies, which minimizes systematics 
when comparing the two epochs. When one uses the same method, it is found that the fraction of fully relaxed E/SO 
does not evolve between z~0.6 and z=0. Compared to other catalogues, an evolution of the fraction of E/SO of the 
order of ~10% is found, but this can be attributed to systematics associated to the different classification schemes 
used. These 10% likely correspond to non-relaxed E/SO, as discussed above. 
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Table 2 

Fractions of morphological type at z=0 in the descendants of the on-going z~0.6 major mergers. Such remnants are expected in a range 
of stellar mass ~lO lo -3-lO ll f M . The error reflect s the uncertainty on the conversion between r-ba nd luminosit y and stellar mass (see 

text). Fractions at z=0 are indicated as derived from Nakamura ct al. (2003), Nair & Abraham (2010), and Dclgado- Serrano ct al. (2010). 

Note that the sum of the fractions derived from Nakamura ct al. (2003) is slightly over 100% because they were derived from the MDLFs. 

Fractions at z~0.6 were determined by [Dclgado-Scrrano ct al. (201(3). 



Morph. type 


Fraction 


Fraction 


Fraction 


Fraction 




at z=0 


at z=0 


at z=0 


at z~0.6 




from N03 


from N10 


from DS10 


from DS10 


E/S0 


28±3% 


31% 


18±5% 


17±3% 


All Sp 










SO/a-Sd 


74±3% 


68% 


72±8% 


31±7% 


Early-type Sp 










SO/a-Sb 


54±4% 


45% 






Late-type Sp 










Sbc-Sd 


20±1% 


23% 






Im / Pec 


2±1% 


<1% 


10±3% 


52±9% 
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